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摘 E 本 文 研究 趋势 项 含有 变 点 且 新 息 为 方差 无 穷 厚 尾 过 程 的 序列 单位 根 检验 问题 ， 通 过 构造 DF 型 检 
验 ， 得 到 了 其 渐 近 分 布 。 为 避免 估计 统计 量 渐 近 分 布 中 的 尾 指数 ， 构 造 subsampling 抽样 方法 来 
确定 统计 量 渐 近 分 布 的 百 分 位 数 ， 同 时 论证 了 subsampling 抽样 方法 的 一 致 性 。 最 后 ， 用 Monte 
Carlo 模拟 证 实 本 文 所 提出 统计 量 以 及 subsampling 抽样 方法 的 可 行 性 。 
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含有 结构 变 点 的 单位 根 检验 问题 是 近年 来 计量 经 济 学 与 统计 学 界 研究 热点 。Perronfl3 提出 
含有 结构 变 点 的 单位 根 检 验 问题 ， 并 构造 DF 型 检验 统计 量 。Banerjee 等 人 四 ，Christiano 向 ， 
Zivot 和 Andrews!) 考虑 结构 变 点 为 内 生变 量 时 的 单位 根 检验 问题 。Kim $6 考虑 新 息 中 存在 
方差 变 点 时 的 单位 根 检 验 问 题 。Nunesl7] 用 LM 型 检验 研究 了 含有 结构 变 点 时 的 单位 根 检验 问 
题 。 另 一 方面 ，Guillaumel8] 以 及 Mittnik 和 Rachevig 指出 许多 经 济 和 人 金融 数据 具有 尖峰 厚 尾 
的 特点 ， 用 方差 无 穷 序列 来 描述 更 能 反应 这 些 统计 特性 。Davis 和 Resnick! 给 出 方差 无 穷 
过 程 的 样本 相关 函数 以 及 样本 协 方差 的 收敛 性 。Chan 和 Tran!) 针对 新 息 为 方差 无 穷 厚 尾 过 程 
情形 下 的 单位 根 检验 问题 ， 使 用 DF 方法 (Dickey-Fuller) 检验 ， 结 果 表 明 直 接 使 用 DF 检验 ， 
检验 的 势 会 有 较 大 的 波动 。Phillipst13] 将 这 一 结果 推广 到 弱 相 依 的 方差 无 穷 情 形 。 在 厚 尾 指数 
以 及 Chan 和 Tran? 中 统计 量 渐 近 分 布 未 知 时 ，Agnieszka 和 Kokoszka0d 应 用 subsampling 方 
法 确定 上 述 DF 统 计量 所 需要 的 临界 值 ， 模 拟 结果 表明 subsampling 方法 有 助 于 提高 检验 的 
势 。 
当 新 息 为 方差 无 穷 过程 时 ， 对 于 在 过 程 的 趋势 项 含有 结构 变 点 时 的 单位 根 检 验 问 题 却 鲜 有 
研究 。 本 文 考虑 了 新 息 为 方差 无 穷 厚 尾 过 程 情形 下 ， 含 有 趋势 变 点 时 的 单位 根 检验 问题 ， 得 到 
此 时 DEF 检验 的 渐 近 分 布 。 然 而 统计 量 的 渐 近 分 布 结构 较为 复杂 ， 且 分 布 中 含有 尾 指数 ， 会 对 
检验 的 势 产 生 影响 。 为 此 ， 本 文 应 用 subsampling 抽样 方法 来 确定 检验 所 需要 的 临界 值 ， 同 时 
可 以 避免 估计 厚 尾 指数 ， 并 证 明了 subsampling 抽样 方法 的 一 致 性 。 数 值 模拟 表明 当 样 本 容量 
较 大 时 ， 其 经 验 势 函数 接近 1， 这 表明 subsampling 方法 是 有 效 的 。 
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2 ”模型 与 假设 
考虑 模型 
1, t < ko, 
Ye = p + ODU; + Bt + pyz_1 + Ut, u=] t=1,---,n, (1) 
0, t>ko, 
其 中 新 息 过 程 {u} 满足 如 下 条 件 : 


假设 1 {wi} 为 独立 同 分 布 的 稳定 分 布 序列 ，Ew: = 0 且 其 特征 函数 为 
exp{—t2/2},，« = 2， 
exp{—|t|"}, 1<<< 2. 


E exp{itui} = | 


检验 原 假设 为 Ho : p = 1， 备 择 假设 为 Hi : |p| < 1。 

对 于 上 述 单 位 根 检验 ， 由 模型 (1) 得 到 p 的 最 小 二 乘 估 计量 p,， 从 而 构造 DF 检验 TT = 
n(pn 一 1)。 由 于 最 小 二 乘 估 计量 pn 与 模型 参数 Jj.、9、B 无 关 ， 因 此 在 本 文 证 明 中 假设 : 

假设 2 yn ,yn 的 生成 过 程 为 y = yi1 +u, t=1, ,n， 其 中 以: 满足 假设 1。 

注 1 为 推导 的 简化 ， 本 节 仅 考虑 新 息 为 独立 同 分 布 的 情形 。 事 实 上 ， 本 文 结论 对 于 更 一 般 
的 弱 相 依 厚 尾 序列 也 是 成 立 的 : w = d(L)er = Do dyat-;， 其 中 系数 do = 1 d(1) #0, Hii 
Eo tldel® <0, 0<6<1Ak， 误 差 {et} 满 足 假设 1。 


3 主要 结果 


在 假设 1、2 成 立时 ， 对 DE 检 验 = n(1 一 pn)， 有 如 下 结果 成 立 。 
引 理 1 车 {yr} 为 单位 根 过 程 ， 即 y = 十 WW， 其 中 多 满足 假设 1 时 ， 则 有 如 下 结论 成 立 


n n 1 
ap Yiu SUN); ag? yom sovi; 


t=1 


n 1 n 1 
n-la-! SO tur 4, ua) -f U(r)dr; n'a; Ymi #, U(r)dr; 
t=1 0 t=1 0 
z eof z a f 
n-*q7! >》 tues 4f rU(r)dr; naz? Y yi sf U?(r)dr; 

t=1 0 t=1 0 


[nà] a M [na] ap 
na,! Yom sf U(r)dr; nlan? Y yi sf U?(r)dr; 
t=1 0 t=1 0 


{nd} 入 [nA] j 入 
n—*a7 So tu- x J rU(r)dr; n`ta7! Ytu — AU (à) 一 f U(r)dr. 
t=1 0 t=1 9 


证 明 在 假设 1 下 ,对 0 <r <1, Hata, S U(r)， 结 合 文献 15] 的 证 明 ， 可 以 证 
明 引 理 1 的 结论 成 立 。 
定理 1 在 假设 1 中 ， 当 原 假设 Ho 成 立 且 变 点 时 刻 i = [n] BAIN, p = 1 的 最 小 二 乘 估 
Th Ôn 满足 
2 a (dc 一 b?)D, — cD, + bD3 + dD, + bDs 
ae (de — bD; —dD; chs + 26D)’ (2) 
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其 中 d=1/12, b= 一 上 ,c= (1-A). 
Dy = 3[u%(a) ~ vay] - va) I | U(r)dr, Ds = (UQ) - AU(1)) I i (r ~ 5) U(ryar, 
= (tua) -f U(r)dr) f (r-5U(r))dr, De= J: U? (rar = (/ 5(ndr) ， 


f Urjar ~ Í i U(r)ar) (300) - f ; U(r)ar), Dr=( f T jarai f i Djdr) 


= ( fe (r)dr — A f ! Ur)ar) (UA) — AUD)), Ds = ( f : (r- 5 )U(rjar)’, 


D2 
Da 
Da 


D = 人 U(r)dr — af U(r)ar) L (r = 5 )U(r)ar, 


SEP U(r) 为 具有 特征 指数 < 的 稳定 分 布 变量 。 
证 明 为 证 明 的 方便 ， 将 对 假设 1 的 回归 分 为 两 步 : 第 一 步 ， 对 模型 各 项 的 常数 作 回归 得 
到 


Yi = Xit + OX2t + pht-1 + Ut, (3) 
其 中 
ic i< 
Y=H--dou, ha =wai-— dow, 
a darn 
i< n+1 1 Š 
Xu=t--9 t=t- -7> Xa = DU; — à, Zi = 二 >》 yer. 
arr ta 
Ew 


a a Xi, oe 4 Xin, 
Y’ = (ĝos ,Ît-1), U’ = (wu, sUn), Z= ‘ 
Xo1, er) Xan, 


第 二 步 ， 对 模型 (3) 进行 回归 ， 得 到 Pp,,， 并 作 如 下 分 解 
fn ~1 = (Y'Y)Y'E —(¥'Y)"12Z[2'(1 - Py)2) 2(1 - Py)B, (4) 


其 中 Py =Y(Y'Y)-1Y'。 若 令 


则 (4) 式 可 以 表示 为 


`, _ datn —b2)Y'E — cn HK + bp JK — dn JL + bn HL 
Pa — (den — b2)Y'Y — dnd? — dn J? — p H? F 2p AJ | 


》 
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结合 dn, Cn, bn, H, J,，K, 工 定义 ， 由 引 理 1 可 以 得 到 (4) 式 右 端 各 项 的 极限 如 下 


nd, =n73(n3/12 —n/12) >d, n?o, = —n~?(1 — A)An?/2 > b, 


Ts 
no, =n "(L-A)An—e, ag'L= y (u — a) $ U(A) MY). 


t=1 


ii ea 1 
az Y'E =a;? Dy Yi)u È z [v2(1) -TYG] vaf U(r)dr 


t=1 


n 1 1 2 
nta; Y'Y = naz? Y (uei — Yu)? 4 | U?(r)dr — ( f U(r)dr) 
t=1 0 0 


n 1 1 
no =n ?a Stu Sale J rU(r)dr 一 : [ U(r)dr 
0 2 Jo 


t=1 


naz = nant Y (t- HA) $ iya) -f U(r)dr 
2 2 2 0 


Tp 入 1 
n-lo-!lJ =n la}? Dy 一 六 1) / U(r)dr — af U(r)dr. 
= 0 0 


对 DF 型 统计 量 ”(pn —1), RANA 


ae (dncn — b2)Y'E — HK +b JK —d,JL + bnHL 
m j (dncn 一 82) — dyJ? — dnJ? — cen H? + 2bn HI 


(de — b?)D, — cD2 + bD3 + dDa + bDs 
(dc EE b2) De = dD7r cDg + 2bDo 


4> 


? 


定理 证 毕 。 

由 定理 1 可 知 ， 当 新 息 过 程 w 为 方差 无 穷 过 程 时 ，DF 型 检验 统计 量 n(p,, — 1) 的 渐 近 分 布 
为 Lévy 过 程 的 泛 函 ， 因 此 临界 值 的 选取 会 依赖 于 对 新 息 过 程 中 的 厚 尾 指 数 ; 的 精确 估计 。 为 
避免 估计 kw， 我 们 用 subsampling 方法 来 逼近 DER 检验 的 渐 近 分 布 ， 从 而 将 抽样 的 经 验 分 布 分 位 
数 作为 DF 型 检验 的 临界 值 ， 具 体 步 骤 为 : 

1) HARE i, = Yk — PnYk—1, k= 2,--+ ,ns 并 将 其 中 心 化 ， 且 记 为 


ue = tp =a Ue 
2) ”基于 上 述 中 心 化 后 的 残 差 过 程 wW， 构 造 n 一 5 个子 过 程 ， 当 = 2,… ,n 一 5 十 1， 第 个 
子 过 程 定 义 为 yo(k) = 0, yi(k) = ub, yo(k) =y (k) = ud yy, ++, yo = Yor + Ukyo 


3) WAP, AAT FILE yolk), --- ,yo(k) 在 模型 (1) 中 回归 所 得 关于 2 的 估计 ， 记 了 2, 大 一 
2,… ,mn 一 b 十 1， 所 生成 经 验 分 布 的 a 分 位 数 为 q9 (a)> A (fn — 1) < ah g(x)» WIE BE 
设 。 
AE n( 一 1) EDA RBA G(r), WIFE subsampling DIER EIN 5 Ai K 
MG = re HT, < z} 有 如 下 的 定理 成 立 。 


定理 2 在 原 假设 下 ， 当 n -00, 5 一 00 且 b/n 一 0 时 ， 对 任意 zx， 均 有 G(z) 5 G(z)。 
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定理 3 在 备 则 假设 下 ， 当 n 一 co, 8 一 oo 且 b/n 一 0 时 ， 对 任意 z， 均 有 
G(x) È H(z)=1. 


注 2 定理 2 表明 在 原 假设 成 立时 ， 用 subsampling 方法 确定 的 临界 值 趋 于 真实 分 布 的 临界 
值 G。 定 理 3 则 表明 在 备 则 假设 成 立时 ， 用 subsampling 抽样 方法 确定 的 临界 值 进行 检验 ， 其 
经 验 势 函数 接近 于 1。 

为 方便 定理 2 的 证 明 ， 先 证 明 如 下 的 引 理 2 至 引 理 5。 

引 理 2 设 随机 变量 序列 Sn, Urn, 1Sk Sm oe > 0, Urm £ Urn Bin = 
Op(1)， 则 对 任意 e > 0, lim sup max P(|Sn Uk,n| > €) = 

引 理 2 的 证 明 参见 文献 1 。 

引 理 3 在 定理 2 的 条 件 下 ， 对 于 上 =2,---,n, u=urtre(n), Hre(n) = Op(nan)- 

证 明 re(n) 的 定义 可 知 

0 


1 n 
rm = Uk — Ue = a ——— > tit 一 Uk 
t=2 


[yk — Ônyr-1ı] — de Pn¥t—1] — Uk 


n 


1 7 
[yk — yk-1 — (Pn ~ 1) yea] — 二 2 (ye — ye—-1 — (An — 1)yt-1] — uk 


1 nr 
pn — 1) (vei Paes Sowa} 
al 


由 于 


= Op(n~"an), Yk-1 = Op(n“*an), 


1 
一 1 


Ta I Dy Yet = OP(n an), n(ô— 1) = Op(1), 
t=2 


从 而 m = Op(n-lan)， 引 理 证 毕 。 
由 于 在 后 续 引 理 的 证 明 中 ， 将 作 类 似 定 理 1 证 明 中 的 分 解 ， 故 为 方便 ， 记 app = 好 /12 一 
5/12，b = —(1 -入 Xp2/2，cok = (1 — A)AB. 


k+b—-1 k+b-1 ETET 1 k+b—1 


Y' Fur = 5 Yi—1Ui 一 De u(i >》, vj); Ak = pD ib- > yj}, 
i=k i=k j=k i=k j=k 
k+b—1 bay Best k+db—1 1 Pes} 
Koy = > (iui - = >», u); Job = >D [pis — (5 D y-1)|， 
i=k j=k i=k j=k 
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k+b-1 k+b-1 k+b— 
Y'El, = y Yi— 1(k)up Jp D Ui Jf Sa Yj- -1(k)], 
i=k i=k jak 
k+b-1 jt 1 k+b-1 k+6-1 
Hik = Xi i[yi- 1(k) 一 (; 2 Yj- i(k) J], Kýr = y (ig -H uw), 
i=k i=k t=k 
k+rb—-1 和 1 天 十 入 b 一 1 1 k+b-1 
k= ia 1(k) — 7 P Yj-1 (r)], R= > (a -3 uw), 
i=k i=k i=k 
k+b—1 TA 
YY), = 5 [oa 1(k) 一 了 ; > Vi-1( w|’, 
i=k 


TH = b “a, ?[(ab,kCb,k — b} p )Y' Evk — Cok Ho,k Kb + bo, Jo,kKo,k 
— abk Jok Lobok + bo Hok Lee], 
T, = bay 7 [(ab,kCb,k — b? p) Y Yb k — abk Jo,kJo,k — Cox Ho, Hbk + 2bo,k Hp, kJo,k], 
类 似 可 定义 72 45 TPP, MAT n(n- 1) = TN,/TP., OO- 1) = Tok [Tok 天 二 2 一 


十 1。 在 这 些 条 件 下 ， 可 以 证 明 。 
引 理 4 在 定理 2 的 条 件 下 ， 对 任意 z， 均 有 lim 十 5 k bt! p(T, <x) 一 G(z)。 


证 明 在 证 明 前 ， 为 方便 ， ciel RN, = TPN ~ TN, RP, ae ~— TP, 注意 到 
事件 (TE, < Xx} 与 事件 {TA + RN, <TR + oN} 等 价 ， 从 而 可 得 


P(TP, < £) < PIA < < TD, + 2e) + PR < -e) + P(rRD, > €), (6) 
P(T, < £) > P(A, < 2TH, — 2€) 一 P(RA >€) + P(xRP < —e). (7) 
车 能 证 明 
Dn eri yi P{|RN | > e} =0, Komisono mar Ae P{|RP,,| > e} =0, Ve>0, 
则 由 上 、 下 极限 的 性 质 以 及 引 理 3 可 得 
1 n—b+1 
lim -一 一 DD P(T}, < £) > G(2), 


n=œ n 一 


AKEE. FOE 
lim sup max  P{|Ri,| > €} = 0， 


2<k<n—-b+1 
对 RP, 可 类 似 处 理 。 记 
[(k+b—1)t] 
Sorlt)=a5! S mn， 
i=k 
则 由 引 理 3， 结 合 5/n 一 0, 1<k&<2， 有 
sup ISo k(l = op(Ż. = =) = op(1). (8) 


poze b+1 
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ia RN, = TRY -TN =U-V+X-Y+Z, $PU, V, X, YZ 定义 如 下 
V= b 40; cor HOEKD — btaz? co pHo k Ko 
X = ba; byk IEL? p — b tas bo k Jok Lbk, 
Y =b ta; aok IELE p — b 10; ab kJo,kLo,k, 
Z= ba; ao kH Len —b 4a; as k Ho kLb,k, 
U = b- 4a; (ab, kCb,k — bgp) Y E? p ~ b ap (ab,kCb,k — bt k) Y’ Bo,k 
= bta; [Y EP, — Y' Eok], 
HERR YEP, t Y’ Eo, 的 定义 可 得 
a,” [Y ER; — Y' Exe] 


k+b-1 k+b—1 
= | DD Yi_1(k)u? 一 5 yi-1(k)u ) 
i=k i=k 
k+b—1 1 k+b-1 1 tbs? 
= (a; > wi ) So (1) + 5Spx(1) — ap? > wiri(n) =- 306 >», r2(n), 
i= i=k i=k 
从 而 由 (8) 式 可 得 
k+b-1 
(a; 5 ui ) Sox(1) = op(1), 
i=k 


结合 引 理 3 可 得 ，ap?[Y'E9. 一 YEsxr] = op(1), AKU = op(1)。 类 似 可 以 证 明 V = 
op(1), X =op(1), Y = op(1), Z = op(1)， 综 合 则 有 


li ax  P{|RN = 0. 
im sup , aX, | {IRF > e} =0 


引 理 5 在 定理 2 的 条 件 成 立时 ， 对 任意 z， 均 有 Var[Gs(7)] 一 0. 
证 明 由 Co(z) 的 定义 有 
1 n—b+1 


Var[Ge(z)] = (n—b)? 5 Var[I{T2 < x}] 
k 


=1 


n—-bn—b-j 


Ces 2_ 2 Cov(HTa < th Tats sz) 


j=l k=1 
= Sı + So. 
由 引 理 4 的 证 明 可 知 ， Si — 0. 将 92 作 如 下 分 解 


b n~b-j 


2 
S2 = wo 和 Cov( TI < £}, {Tkt < z}) 


n—b n—b—j 


tara LL Corll < ah Bers <3) 


j=b+1 k=1 


= Sor + Soe. 
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由 b/n 一 0 可 得 S21 一 0， 下 证 Soo — 0. 由 引 理 4 的 证 明 过 程 可 知 


P{IRE| > €} + PileRp,| > ef] — 0 


= max 
ý Ba | 


EX pt = PT < zT +2e), p; = P(T, < xT, — 28), 
Pij = slits aT, + 2€, Trea S ITP 4; + 2€), 
Pag = P( TN. < 2TH, + 2€, TS 2TH. — 2€). 
由 (6) 45 (7) 式 可 得 
Pn — ôn < P(Tòk < 2) < prt + Sn, 


Pnij 一 Ôn < P(TD, < T, TÈ ,天 十 了 > < z) < Phy + 26n, 


因此 
Cov( {Tk < £}, {Tp na; <& 7}) < pt; (Pn)? + 26n + Oa, 
从 而 
2 S 
S22 < Goa? p2 (n-b-j- 1){p}; — (py)? + 25n + 62}, 
j=b+1 
故 


+)2 -)2 +_— y + ~_} a 
lim sup S22 < 2[ (p 7- (p7) ], p*= lim ph, p` = lim p. 


H e > 0 的 任意 性 可 得 lim sup S22 = 0. 


定理 2 的 证 明 ”为 证 Gs(z) => G(xz)， 仅 需 证 明 E[Gs(z) 一 G(x)] — 0, Var[G,(x) - G(z)] 一 
0。 而 由 引 理 4、5 WLS EÂ (x) - G(x)] 一 0，Var[Gs(z) - G(z)] 一 0， 从 而 定理 证 毕 。 

为 证 明定 理 3， 先 证 明 如 下 3 个 引 理 。 

引 理 6 在 定理 3 的 条 件 下 ， 有 wh = wx 十 Yk 十 7(n)， 此 处 


T n 
1 k 
Yk = (bn — P)yr- nn) = -~—— 2- (Pn =p) De } 
t 一 2 t=2 


Ai. n(n) = Op(n an)。 
证 明 由 定义 


u? 一 uz = ù- 3 tik 一 Uk 
n 


Á 1 n 
= [yk — ÔnYe-1] — = 2 luk — Pye] — ur 
k=2 
n 


% 1 A 
[ur — pyk-1 — (Ôn ~ P)ur-1] — —— 2 [Yk — pyr-1 + (P — Pn)] — ur 
k=2 


— (Ôn — P)Yk-1 — —{ Du- (in - ») ues} 
k=2 


= yk + 7(n), 
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而 由 引 理 1 可 知 ，n(n) = Op(n-1an)， 从 而 引 理 证 毕 。 
引 理 7 在 定理 3 的 条 件 下 ， 对 任意 z， 均 有 


n—b+1 


; 5O P(E, <2) > H(z). 
k=2 


证 明 与 引 理 4 的 证 明 类 似 ， 仅 需 在 备 则 假设 下 证 明 下 面 两 式 成 立 


lim 
n=% n 一 


N = = U. 
limsup ,max oe P{ |R| > € } =0, limsup max, P{|IRM|>e}=0 


以 下 仅 证 明 第 一 项 ， 第 三 项 可 类 似 证 明 。 在 引 理 4 与 引 理 6 的 相关 记号 下 ， 有 
RY, =T -TA =U-V+X-Y +Z, 
其 中 
U = bta; (ab, kCo,k ~ bgp) Y¥ Ep, — 0 as (ab, kCo,k — 05,4) ¥ Eo,k 
= O(b-*)a,7[Y’ Ep, —Y’ Bk], 
a, [YER p — Y' Eaa] 


k+b—1 k+b-1 


= on 人 ( Ð uat- Y wali) 
i=k i=k 
k+b—1 k+b—1 kto—1 
(a Ee E wno E o] 
i=k i=k 
k+b~1 k+b~-1 


+a,” 2 ui( (n(n) +) ~ 505? 2 (n(n) +%) 


k+b—1 k+b—1 k+b-1 Kaia 1 k+b-1 


= {(a;" DS us) (a5? Yan m) +5 =a,” SS nf n) +a,” > uin(n) — 505? 2 n(n) } 
i=k i=k i=k i=k 
k+b—1 k+b—1 
+{ — a; D nlai ~ 595° DD #} 
= = 

k+b—1 1 k+b-1 k+b-1 Rom 1 A 1 

+{ a5? > nn) + 505" D of - (43? > ui Jap’ 3 yi ta,” D uy } 
i=k i= i= i=k 


UR + UR? + UM 
由 于 当 n 一 off, b/n 一 co， 结合 引 理 6 可 得 


i 


N1 
=0, 
P fra a [Uk (>e) 


而 由 文献 上 4 中 引 理 6 的 证 明 可 得 


N2 p 
limsup | _ max a Use| >e} =9, limsup max, P{ Ub l>e} =O, 
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Re ER, WAU = op(1)， 类 似 可 证 V = op(1), X = op(i), Y = op(1), Z = op(1)， 从 而 


lim sup P aE P{|RM| >e}=0, 
故 引 理 证 毕 。 
引 理 8 在 定理 3 的 条 件 下 ， 对 任意 z， 均 有 Var[Go(z)] 一 0。 
证 明 证明 与 引 理 5 相同 ， 故 略 。 
定理 3 的 证 明 为 证 Gy(z) > G(z)， 仅 需 证 明 E[Gs(z) - G(z)] 0, VarlGo(z) — G(z)] 一 
0。 而 由 引 理 7、8 可 以 得 到 [Gs(z) 一 G(z)] +0, Var[G,(x) - G(z)] 一 0， 从 而 定理 证 毕 。 


4 ”模拟 研究 
本 节 用 模拟 研究 来 验证 本 文中 subsampling 方 法 的 有 效 性 。 考 虑 如 下 数据 生成 过 程 (DGP) 
1, t< ko, 
ye = p + ODU; + pyz—1 + Ue, om -1 t 一 1 (9) 
0, t> ko, 


yo = 0， 其 中 参数 jy = 0.1, 0=0.2, DU, 如 前 所 定义 ， 为 厚 尾 指数 上 = 1.14 的 独立 随机 变量 
序列 。 

表 ]1 为 在 原 假 设 下 ， 即 p = 1， 样 本 量 n = 1000 时 ， 直 接 模拟 重复 5000 次 试验 所 得 百 分 位 
数 ， 表 2 为 样本 量 n = 1000 时 ， 子 样本 量 为 m = [lIn(n)?]， 重 复 抽 样 5000 次 时 所 得 百 分 位 数 。 
对 比 表 1 与 表 2，subsampling 方 法 确定 的 百 分 位 数 波动 较 小 ， 且 计算 运行 时 间 较 短 。 表 3 与 
表 4 分 别 为 在 备 择 假设 下 ，p = 0.7， 利 用 表 1 与 表 2 中 的 分 位 数 计算 所 得 经 验 势 函数 值 ， 重 复 
试验 次 数 为 5000 次 。 对 比 表 3 与 表 4，subsampling 抽样 所 得 临界 值 的 势 与 直接 模拟 所 得 临界 
值 相近 ， 同 时 ， 当 样本 量 n 逐 渐 增 大 时 ， 两 种 方法 的 经 验 势 函数 均 接 近 于 1。 


表 1: 当 k 二 1.14 时 ，Monte Carlo 模拟 所 得 渐 近 分 布 (2) 的 百 分 位 数 


1% 2.5% 5% 10% 90% 95% 97.5% 99% 
A=0.1 -49.6916 -33.0439 -25.8829 -20.6187 -4.8297 -3.7336 -2.9012 -1.8471 
0.2 -55.8427 -36.2811 -26.9334 -21.0022 -5.2803 -4.0516 -3.0127 -1.5167 
0.3 -57.3373 -36.4277 -26.9860 -21.1285 -4.8980 -3.5957 -2.4422 -0.8428 
0.4 -53.4348 -36.7738 -27.8241 -20.9102 -4.3967 -3.1490 -2.1149 -0.6182 
0.5 -54.7308 -35.7396 -26.9451 -20.8586 -3.9839 -2.9013 -1.7961 -0.4268 
0.6 -52.9663 -38.3679 -28.2609 -21.2167 -4.4016 -3.1920 -1.9892 -0.6271 
0.7 -64.5498 -37.2638 -28.0128 -21.5163 -5.0157 -3.8594 -2.6479 -1.0842 
0.8 -60.3578 -36.6004 -27.8867 -21.4129 -5.3179 -4.0841 -3.0800 -1.8164 
0.9 -52.3821 -34.7184 -26.6580 -20.3535 -4.8968 -3.7787 -2.8314 -1.4228 
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表 2: 当 k& = 二 1.14 时 ，subsampling 方 法 模拟 所 得 渐 近 分 布 (2) 的 百 分 位 数 
1% 2.5% 5% 10% 90% 95% 97.5% 99% 
A=0.1 -45.1578 -25.7280 -21.4191 -16.7911 -5.0268 -3.6300 -2.3096 -1.7618 
0.2 -42.2333 -29.1284 -23.0033 -18.2148 -5.4403 -4.3610 -3.4893 -1.3655 
0.3 -43.3304 -30.2418 -23.0136 -18.6540 -5.2126 -4.2637 -3.6270 -2.1643 
0.4 -44.3758 -30.7756 -24.0014 -18.8300 -4.6412 -3.8584 -2.6047 -1.6014 
0.5 -45.1691 -36.9219 -30.4722 -22.9231 -4.7443 -4.0289 -2.9999 -1.7521 
0.6 -45.1691 -36.9219 -30.5058 -22.9231 -4.7055 -3.9355 -2.7280 -1.0672 
0.7 -45.1691 -36.9219 -30.4722 -22.5792 -4.8835 -3.8678 -2.7280 -1.0672 
0.8 -45.1691 -32.2741 -28.8533 -21.5030 -4.8835 -4.0178 -2.6694 -1.5482 
0.9 -45.1691 -31.6323 -26.5962 -20.7075 -4.8404 -3.8822 -2.5633 -0.9671 
3: 当 乒 =1.14 时 ， 不 同样 本 容量 时 Monte Carlo 模 拟 所 得 经 验 势 函数 
n = 100 n = 200 n = 300 
入 1% 2.5% 5% 1% 2.5% 5% 1% 2.5% 5% 
0.1 0.0536 0.6524 0.9516 0.9764 0.9998 1.0000 0.9996 0.9998 0.9998 
0.2 0.0126 0.4602 0.9430 0.9162 0.9996 90.9998 0.9994 0.9998 0.9998 
0.3 0.0092 0.4474 0.9394 0.8970 0.9994 0.9996 0.9992 1.0000 1.0000 
0.4 0.0186 0.4048 0.9272 0.9526 0.9990 1.0000 0.9996 1.0000 1.0000 
0.5 0.0156 0.4782 0.9364 0.9364 0.9994 0.9998 0.9996 1.0000 1.0000 
0.6 0.0224 0.3060 0.9182 0.9564 0.9980 0.9994 0.9996 1.0000 1.0000 
0.7 0.0044 0.3866 0.9204 0.5674 0.9998 0.9998 0.9970 0.9998 0.9998 
0.8 0.0058 0.4246 0.9180 0.8024 0.9998 1.0000 0.9986 1.0000 1.0000 
0.9 0.0392 0.5296 0.9380 0.9598 0.9990 0.9998 0.9998 1.0000 1.0000 
表 4: 当 一 1.14 时 ,不 同样 本 容量 时 subsampling 方 法 所 得 经 验 势 函数 
n = 100 n = 200 n = 300 
入 1% 2.5% 5% 1% 2.5% 5% 1% 2.5% 5% 
0.1 0.1882 0.9544 0.9902 0.9968 1.0000 1.0000 1.0000 1.0000 1.0000 
0.2 0.1542 0.8884 0.9864 0.9970 1.0000 1.0000 1.0000 1.0000 1.0000 
0.3 0.1280 0.8618 0.9868 0.9936 0.9994 0.9996 1.0000 1.0000 1.0000 
0.4 0.0932 0.8300 0.9792 0.9938 0.9994 0.9996 1.0000 1.0000 1.0000 
0.5 0.0838 0.3876 0.8464 0.9922 0.9990 0.9998 0.9998 0.9998 0.9998 
06 0.0818 0.3974 0.8388 0.9920 0.9992 0.9996 0.9998 0.9998 0.9998 
0.7 0.0886 0.4174 0.8500 0.9924 0.9996 0.9998 1.0000 1.0000 1.0000 
0.8 0.0896 0.7396 0.9000 0.9914 0.9994 0.9996 1.0000 1.0000 1.0000 
0.9 0.1050 0.7596 0.9430 0.9914 0.9998 0.9998 1.0000 1.0000 1.0000 
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Subsampling Procedures for the Heavy-tailed Unit Root Test with 
Structural Change 
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Abstract: This paper deals with the unit root test for the series with the infinite-variance innovations 
and the change in trend. The DF-type test is extended and the asymptotic distribution is obtained. 
In order to get the percentiles of the DF-type test, subsampling procedures are designed, and the 
consistency of the subsampling procedure is established. Finally, Monte Carlo simulations demonstrate 
the validity of the proposed test and the subsampling procedures in this paper. 
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